Evidence for nonmonotonic magnetic field penetration in a type-I superconductor 
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Polarized neutron reflectometry (PNR) provides evidence that nonlocal electrodynamics governs the 
magnetic field penetration in an extreme 1ow-k superconductor. The sample is an indium film with 
a large elastic mean free path (11 (im) deposited on a silicon oxide wafer. It is shown that PNR 
can resolve the difference between the reflected neutron spin asymmetries predicted by the local and 
nonlocal theories of superconductivity. The experimental data support the nonlocal theory, which 
predicts a nonmonotonic decay of the magnetic field. 

PACS numbers: 74.20.-z, 74.25.Ha, 78.70.Nx 



In this paper we pose and answer experimentally the 
following fundamental questions. Are nonlocal electrody- 
namics effects measurable in superconductors? Can the 
nonmonotonic decay of magnetic field penetration pre- 
dicted by the nonlocal theory be observed? To what ex- 
tent can Polarized Neutron Reflectometry (PNR) resolve 
the difference between local and nonlocal diamagnetic re- 
sponses expected for strongly type-I superconductors? 

Nonlocality is a key concept of superconductivity the- 
ory, but its experimental verification is still not estab- 
lished. In the Meissner state, a magnetic field ap- 
plied parallel to the surface located at z — causes 
the magnetic induction B(z) to penetrate over a depth 
A = -B(O)- 1 / B(z)dz. In the London (local) limit, 
appropriate to most type-II superconductors, B(z) oc 
cxp(— z/Xl), where Xi is the London penetration depth. 
In 1953, to explain the variation of A in type-I supercon- 
ductor Sn caused by adding In, Pippard proposed that 
the current density is related to the average of the vector 
potential over a region of size £o (the Pippard coherence 
length) m • The smaller the Ginzburg-Landau parameter 
k, the more important this nonlocal effect. 

Nonlocal theory predicts that B(z) deviates from a 
simple exponential decay. B(z) is nonmonotonic and, 
moreover, changes sign at a specific depth In the 
pure limit (£o t, where I is the elastic mean free 
path) B(z) is determined by the intrinsic parameters 
Xl (T = 0) and £ , an d by the temperature T. The devia- 
tion is most significant in "extreme" type-I superconduc- 
tors (k < l/\/2), such as Al (« « 0.01) and In (0.06). 
For these the results of the Pippard theory are identi- 
cal to those of the Bardcen-Cooper-Schrieffer theory Q. 
On the other hand, the local approximation is safe for 
> 1.5 Q, i.e., marginal nonlocal effects are predicted 




for some type-II superconductors such as Nb (k w 1). 



50 100 150 200 250 300 350 400 

z (nm) 

FIG. I: Magnetic-induction profiles B(z) in a semi-infinite In 
sample. The dashed (solid) line corresponds to the local (non- 
local) relation between current density and vector potential. 



B(z) in In, calculated in local and nonlocal approaches, 
in the pure limit with £o = 0.38 and Xl(0) = 0.025 
/xm for T = 1.8 K, is shown in Fig. 1. Details about 
the formalism can be found in Ref. 3. In the nonlocal 
approach B(2Xl) / B(Q) w 1/e. The sign reversal is ex- 
pected at z sa 5.5 A/,, and the amplitude of the reversed 
field is at most about 3% of the field at the surface. 

An observation of sign reversal was reported in Ref. 
5. An external AC magnetic field H with amplitude up 
to 30 Oe was applied parallel to a hollow cylindrical Sn 
film about 2 /im thick, and a strongly attenuated (10 s 
times) signal with reversed phase was detected inside the 
cylinder at T = 2.88 K and H w 25 Oe. This phase 
difference was interpreted as a change of sign in the pen- 
etrating field. However, this is questionable because the 
phase difference drops back to zero at a larger (30 Oe) 
field, whereas the critical field H c at 2.9 K is 115 Oe Ira]. 
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Nowadays B(z) can be measured directly using po- 
larized neutron reflectometry (PNR) 7] and low-energy 
muon spin rotation (LE-//SR) [H techniques. We com- 
ment briefly on the latter before focusing on the former. 

In the LE-/iSR technique positive muons polarized per- 
pendicularly to the applied field are implanted in a sam- 
ple over a distance determined by the muon energy. B(z) 
is obtained from measuring the Larmor precession fre- 
quency of the muon spins at stopping distance. In other 
words, the implanted muons serve as tiny sensors of the 
magnetic field inside the sample. However, in practice the 
muon precession is progressively damped with depth due 
to a broad distribution of stopping distances [§]. From 
our prospectus, this is the main difficulty in applying the 
LE-/iSR technique to fields with a sharp profile. 

Recently, the LE-/iSR technique was used to measure 
B(z) in Pb, Nb, and Ta [9]. Most interesting is the obser- 
vation of a non-exponential shape of B{z) for all studied 
metals. The nonlinearity of the semi-log plots for B(z) 
is marginal, which is exactly what should be expected 
theoretically in view of the fairly high k of the studied 
samples. For example, k of pure Nb (residual resistivity 
ratio RRR = 1600) is 1.3 at 3 K and 1.0 at 7 K [10]. How- 
ever, in Ref.[^ n of less pure Nb (RRR=133) is reported 
to be 0.7 at 2.96 K and 7.6 K. This and some other in- 
consistencies with well established literature data suggest 
that the muon probing results may contain some hidden 
uncertainties. Therefore, additional experiments would 
be worthwhile, in particular on low-n superconductors. 

The PNR technique is based on the change of the neu- 
tron index of refraction in a magnetized medium. When 
a collimated neutron beam polarized along the magnetic 
field is incident on a flat, laterally uniform sample under 
a grazing angle, its specular reflectivity R is determined 
by the profile of the neutron scattering potential below 
the surface. R is measured versus momentum transfer 
Q = 47r sin 9/ X n , where 9 is the angle of incidence and 
X n the neutron wavelength. The scattering potential con- 
sists of a nuclear and a magnetic part, which results in 
different reflectivities R + and R~ for neutrons with spins 
polarized parallel (up) and anti-parallel (down) to the 
applied field, respectively. Direct information about the 
sample magnetization is obtained by combining R + and 
R~; the combination s = (R + — R~)/(R + + R~) is the 
spin asymmetry. B(z) can be found by fitting s(Q) data 
with s(Q) calculations based on theoretical models for 
B{z). PNR has been applied for measuring the penetra- 
tion depth and for detecting surface superconductivity in 
Nb dUd, high-T c cuprates [H Q El and Pb PS 03 . 

The nonlocal effect in B{z) measured with PNR was 
discussed in Ref s . HT1. fl2l fl6l and li~7l Although some devi- 
ation from exponential decay was noticed in Rcfs. and 
liH no solid confirmation of the nonlocal theory was ob- 
tained. The authors of Ref.[l3 correctly pointed out that 
experiments with 1ow-k type-I superconductors are de- 
sirable to verify nonlocality, but their overall conclusion 
was that PNR is incapable of detecting nonlocality in 
any superconductor. In this light it is very interesting to 



reassess the problem of nonlocality with state-of-the-art 
PNR applied to a 1ow-k material such as In, since during 
the last decade many neutron source facilities have sig- 
nificantly progressed in neutron flux, neutron optics and 
detector technology. 

The design of the sample for the PNR study is based 
on the following requirements. The irradiated surface 
must be flat and possess minimal possible roughness. The 
sample must be thick enough to have the same electro- 
magnetic properties as the bulk material. Degradation 
of the surface quality with increasing thickness limits the 
film thickness. Neutrons reflected back from the sub- 
strate should have a negligible effect on the reflectivity 
in a region close to the critical edge of total reflection, Q c , 
where the reflectivity is most sensitive to the magnetic 
properties. 

Two approaches can meet these requirements. One 
is to deposit a thick film on a flat substrate that reflects 
least. This can be achieved if the neutron refraction index 
of the substrate is larger than that of the sample. This 
approach was taken in the experiments on Nb [lll.ll2j] and 
Pb 0, In fact) this was the only option, in view of 
the negligibly small absorption of neutrons in Nb and Pb. 
However, In is a strong absorber, which enables one to 
rely on substrates with a refractive index smaller than 
that of In, provided the thickness of the indium film is 
properly optimized. In this approach a second plateau or 
"hill", associated with total reflection from the sample- 
substrate interface, is expected in the reflectivity curve 
R(Q). This should yield additional information about 
the sample structure. Modelling shows that an indium 
thickness of 2.5 /xm is appropriate. Such a sample was 
fabricated in the present work. 

High purity indium (99.9999%) was deposited by ther- 
mal evaporation on the polished side of a silicon oxide 
wafer at room temperature. The substrate size was 2x2 
cm 2 x 1 mm. The base pressure and the evaporation rate 
were 4xl0 -8 mbar and 60-70 A/s, respectively. The 
nominal film thickness, as recorded by a quartz monitor, 
was 2.5 /im. Several smaller area samples were simulta- 
neously fabricated for the film characterization. 

The root-mean-square (rms) surface roughness a 
probed with an atomic force microscope (AFM) yielded 
2.0, 6.7 and 8.0 nm at the scale of 1, 5 and 10 /im, re- 
spectively. A scan range up to 10 /xm was not sufficient 
to reach saturation of the roughness. Consequently, 8.0 
nm is a lower bound on the roughness at the scale of 
the neutron coherence length (ss 100 /im [18|). In our 
simulations, effects due to surface roughness are mod- 
elled using Nevot-Croce factors p}, where the roughness 
is characterized by <r; it was allowed to vary to fit the 
experimental data. 

Another parameter associated with the sample surface 
is the thickness of the indium oxide film. When exposed 
to air, In, like its neighbors in the Periodic Table, Al and 
Ga, instantly forms a protective oxide layer. A surface of 
indium in air remains lustrous for years. This suggests 
that the oxide layer is very thin, perhaps of the order 



3 



of a few monolayers, and should not affect the neutron 
reflectivity. This is consistent with the negative result of 
Rutherford backscattering measurements performed on 
our sample: no oxide film has been detected. 

The electromagnetic properties of the sample were 
characterized by measurements of the DC magnetiza- 
tion M and the electrical resistivity. The shape of the 
M(/f)-curves is typical for type-I superconductors. The 
obtained phase diagram H C (T) agrees well with the lit- 
erature data T c of our sample (3.415 K) matches the 
tabulated value of 3.4145 K [ijf , and RRR = 540. Corre- 
spondingly, i w 11 fjxa is much larger than £n- Therefore, 
our sample is a type-I superconductor in the pure limit. 

PNR experiments were performed both on the RE- 
MUR reflectometer [20| at the Joint Institute for Nuclear 
Research (Dubna) and on the CRISP instrument [21| at 
ISIS (Oxford). Both sets of measurements confirm that 
splitting of the R + (Q) and R~(Q) curves is achievable for 
our sample. The ISIS data, which are the most detailed, 
allow a quantitative analysis to which we now turn. 

CRISP operates with a spin-polarized polychromatic 
pulsed neutron beam. The angle of incidence and the 
instrumental resolution AQ/Q were set to 0.24 degrees 
and 3%, respectively. 

The reflectivity in the Meissner state was measured at 
T = 1.8 K and at magnetic fields of 77, 140, 166, and 194 
Oe (H c (1.8 K) = 205 Oe). The obtained data sets are 
shown in Fig. 2. The i?(Q)-dependencies exhibit a hill 
caused by the total reflection from the substrate. The 
splitting between R + and R~ is clearly visible near Q c ; 
different magnitudes of the error bars are due to different 
times of exposure. The data obtained at 77 and 166 Oe 
have the smallest statistical error and will be used for 
further discussion. 

The data for the reflectivity in the normal state are 
shown in Fig. 3. Solid curves are simulations, in which 
the sample was represented by a pure In film on a Si02 
substrate. In the simulations the angular beam resolution 
was allowed to vary due to the unknown uncertainties of 
the instrumental resolution and of the geometrical factor 
(as only part of the beam covers the sample). 

The simulation curve near Q c is mostly controlled 
by the resolution (see also [l6j]). The next segment, 
down to the foothill, is determined by the roughness of 
the sample surface. The location of the ascending part 
(0.011 < ^(A- 1 ) < 0.014) is governed by the film thick- 
ness. The curve segment following the hill is determined 
by the substrate scattering properties. No attempts were 
made to achieve a better fit for that segment, because 
there the spin asymmetry is indistinguishable from zero. 

The best fit (Fig. 3) was obtained for the model sample 
with a — 14 nm and AQ/Q — 2.5 %. Fitting the ascend- 
ing part enables one to determine the film thickness in 
situ. The statistical error of the reflectivity data in this 
region being ± 5%, the thickness was found to be 2400 
± 30 nm, in agreement with the nominal thickness of 2.5 
/im. These parameters were further used for simulating 
the spin asymmetry. Attempts to introduce an indium 
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FIG. 2: Reflectivity of polarized neutrons measured in the 
Meissner state. Q c is the momentum transfer for total neu- 
tron reflection from the outer surface. The scale is shown for 
the data at H = 194 Oe; the other data have been shifted for 
clarity. 
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FIG. 3: Neutron reflectivity at T = 4.6 K. Curves 1 and 
2 are simulations for a film thickness of 2.40 and 2.50 /im, 
respectively. The inset shows the data for the full range of 
Q-vafues. 



oxide layer on top of the sample yielded no reasonable fit 
for any appreciable thickness (> lnm) of the oxide layer. 
This is consistent with our expectation that the indium 
oxide layer does not affect the neutron reflectivity. 

For our simulations of the reflectivity in the Meissner 
state, the magnetic field profiles shown in Fig. 1 were as- 
sumed on both sides of the sample. The spin asymmetry 
data for fields 77 Oe and 166 Oe, along with simulations 
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FIG. 4: Spin asymmetry at T = 1.8 K and H = 77 Oo (a) 
and 166 Oe (b). The curves are simulations performed within 
the local (dashed line) and nonlocal (solid line) approaches. 



for the local and nonlocal field distributions, are shown 
in Fig. 4. 

For field 77 Oe (Fig. 4a), the results of the "nonlocal" 
simulation fit the experimental data somewhat better, 
but no clear discrimination between the local and nonlo- 



cal approaches is possible due to insufficient accuracy of 
the data at this field. A significantly clearer distinction 
is apparent for field 166 Oe due to the larger amplitude 
of s(Q). As can be seen from Fig. 4b, the quality of the 
fits, with B{z) calculated in the local and nonlocal ap- 
proaches, is different. The nonlocal simulation fits the 
experimental data definitely better, ft is worth stress- 
ing that no adjustable parameters have been used for the 
simulations of spin asymmetry. 

In conclusion, nonlocal electrodynamics effects are 
measurable in extreme type-I superconductors. State-of- 
the-art PNR measurements performed on 1ow-k super- 
conductor In, combined with simulation, unambiguously 
support the nonlocal theory and at the same time demon- 
strate consistency with the literature data for Al(0) and 
£o- Consequently, evidence has been gathered for the 
nonmonotonic decay and sign reversal of the penetrating 
magnetic field predicted by the nonlocal electrodynamics 
approach. 
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